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TECHNICAL MEMORANDUM NO. 1157 

CONICAL TIPS IN SOEERSONIC PLOW* 
By V. Hantzscho and H. Wendt 


In the case of conea In axially symmetric flow of aupersonlc 
velocity, adiabatic compression takes place between shock wave and 
surface of the cone. Interpolation curves between shock polar s 
and the surface are therefore necessary for the complete under- 
standing of this type of flow. They are given In the present 
report by graphical-numerical Integration of the differential 
equation for all cone angleB and airspeeds. 


INTRODUCTION 


The assumption that a velocity potential exists in gas flows 
with supersonic velocity Is compiled with In few specific cases 
only- The expansion of a gas around a corner (reference l) , a 
closed lntegrable example of .a supersonic flow, and the flow 
against a wedge are common knowledge* It Is true that a compressi- 
bility shock occurs on the latter, but the flow after the shock 
Is a parallel flow again, at least In the cases where the shock 
begine at the tip of the wedge- In 1929 Busemann (reference 2} 
discussed the caBe of conical tips in axial flow. In this Instance, 
when the cone angle for fixed Mach number is not too great, the 
flow after the shock Is the simplest case of a more general type, 
the so-called conical fieldB, subsequently discussed In greater 
detail by Busemann (reference 3). These are (not necessarily 
rotation free) dimensionless flows which are reflected on them- 
selves by similar distortion from a certain point. AH straight 
llneB through the center of similarity are carriers of equal 
states of gas (speed, pressure, density) • The axial conical field 
following a parallel . flow occurs, aside for flow against conical 
tips, also for a certain type of compression nozzle. 

The present - report -is Intimately linked to B us smarm 's axial 
conical fields, but deals only with cones. The exact data and 
the further Interpretation of Busemann' a apple curves (Apfelkurven) 
are given. 


*"Mit Oberschallgeschvlndlgkelt angeblasene Kegelspitzen, " 
Jahrbuch 19^*2 der deutschen Luf tf ihrtf orschung, pp. ifiO - 190. 
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THE WEDGE IS SUPERSONIC FLOW 


If a compressibility shock occurs In a two-dimensional 
stationary gas flow, the following equations from the continuity, 
the momentum balance, and the energy law are applicable: 


Pl w ln “ P2"2n 

(1) 

Pi + Pl w ln 2 ■ Vz + 02 v 2n 2 

(2a) 

p l w ln V lt " p 2 w 2n w 2t 

(2b) 

il + 2 « ig + J v 2 2 

(3) 


Subscripts: 

1 quantities before shock 

2 quantities after shock 

Quantities v n and w t are the velocity components normal and 
parallel to the shock surface, p the density, p the pressure, 
and 1 the heat content (fig* l) • Equation (2b), together with 
equation (l), indicates that w^ a w^g, hence, that the tangential 

component of the velocity remains unchanged at the shock. 


The shock process is preferably followed in holographs with 
cartesian coordinates u, o. The states (ug,ug) obtainable' 
by shook from a. certain Initial speed located, say in the 
u-azls (u^), lie on the closed branch of a cartesian sheet of 

the shock polars, (fig. 2) whose equation for an ideal gas 
is readily derived from equations (l) to (3): 



=. (ui - ug 



°P 

a* is the critical velocity, it * r— the ratio of the specific 
■heats at constant pressure and volume* The normal on the length 


( 4 ) 
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also' shows far fixed 


oonnecting the ’state befpre.iand after the shock, a chord, obviously 
Indicates the direction of the shock surface with respect to the 
flow, because at the shock only the* normal^ fe canponent is subjected 
to a discontinuity. Thus the floy against „a wedge vith the angle 
£6 in figure 3 corresponds -to the conditions represented in fig- 
ure £■ The limiting value toward which the angle tends at .decreasing' 
Intensity of compressibility shock is, therefore, the Mach angle- TM 
upper limit far 7 is 9Q° and appears with a straight ccmpreasir 

bility shock (“2 '* ”2 ’ °) - Figure 2 

Mach' number a maximum deflection, that is, a greatest possible 
wedge angle. Far wedges with smaller angle two compressibility , 
shocks are possible, one of which generally leads to superso3U.es 
and one to subBonics* The practioally occurring solution, denoted 
as stable, is the supersonic solution, which, as demonstrated by 
G- Guderley can also be proved experimentally. Since the contact 
point of the tangent from the zero point of the hodograph to the 
shock polar does not, in general, coincide with the Intersection 
point A of the circle of the critical velocity, but still lleB 
in the subsonic zone, there are two subsonic solutions for wedge 
angles on the, admittedly, very small arc between A and B. 


For a wedge in a flow at angle greater than the maximum, a 
curved shock area before the wedge results. The flow after the 
shock is no longer potential flow, the entropy Jump along each 
streamline is different (fig* 4). All points of the shock polar 
are realized. Far from the wedge the compressibility shock runs 
out in a Mach wave} this corresponds to the supersonic end of the 
shock polar s* On the line of symmetry before the wedge 1 b a straight 
compressibility shock, that is, the subsonic end of the polars- 
Between the two states lie in steady sequence the other points, 
particularly point A of the critical velocity and point B of the 
maximum deflection. On the subject of separation of the shock 
area from the tip, L. Crocco (reference 4) established, in an 
exact theoretical analysis of slightly curved wedge surfaces, 
that this separation must take place before reaching the maximum 
angle, which is, to be exact , for a point of the shock polars 
between A and B. One essential condition being, of course, that the' 
curvature at the wedge edge.is other then zero. On considering 
Crocoo's critical point on the shock area separated from the body, 
the streamline there can start only with the. curvature zero* 


Lastly, the variation of the pressure difference pg - p^ 

■ . _ . p 1 

referred to dynamic pressure and the' ratio of 
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the tank pressure before and after the shock follows as 


*g ~- p - 2 






(5) 


/ ' .[•(«« + 1) ban (7 - 8) 

5 (k ¥ l) tan 7 - (K ■ 


- (ft - l) tan 7 ] 


JL 

K-l 


l) tan (7 " 8) 


( 6 ) 

tan ( 7 . - 8} 

J 


with equation (6) tho entropy Jpmp 

As ** R In r (H “ gas constant) 

'Po - 

ls known. . 

« 

. THE CONE IN SUPERSONIC FLOW ‘ 


For a solid of rotation in straight flow, when the shock sur- 
face adheres to the cone tip, the flow after the shock is a conical 
field- The shock surface is a coaxial cone of rotation (fig. 5) • 

The: compressibility shock on the entire cone envelope has the same 
strength. In all points of each straight line through the cone 
origin the velocity is constant in magnitude and direction; the 
state of gas,, the same. The shock .conditions (l) to (2) are 
exactly the seme as for two-dimensional flow, when,- as before, 
w t and w n denote the velocity components parallel and at right 

angles to the shook, area. (Owing to the axial symmetry, w-t and v n 
lie in an axial section of the cone.) Their erection merely requires 
the plane tangential element to the shock areai The potential 
flow following the shock is naturally no longer a parallel flow in 
the axial section it would be incompatible with the continuity 
condition - the streamlines rather become asymptotically narrower 
(fig. 6). Figure 7 shows the conditions in the hodograph for one 
axial section, with u and n denoting the velocity components 
parallel and at right angles to the axis of rotation. Because of 
the conical character of the flow in the hodograph all streamlines 
are represented by the. same curve. The velocity before the shock 
is along the axis of* the abscissa (uj.) • At the shock cone with 
the cone angle 27 tho streamlines are deflected by the angle 8, 
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as seen from the shock polar related to u^> The nev state ug, ^ 

Is transferred by adiabatic compression to the state u^, ^ at the 

cone envelope* Angle p is then half the included angle of the cone* 
A closer examination of the conical field Indicates (compare refer- 
ences 2 and 3) that the normal in a certain point 5 1 of the meridian 
curve runs parallel to the straight line through the cone tip, which 
carries the gas state defined by P’ , as exemplified in figures 6 
and 7 for the state u 1 ., o' with the normal of the slope ft'* Far 
the end point u^, direction of velocity end slope of the normal 

are coincident. Owing to this arrangement of streamline and 
meridian line, the specification of a meridian lino characterises 
a cone in axial flow* To obtain a survey over all conical flows, 
either the cone opening angle at the same aircpeed, hence the shock 
polar s, or, which 1 b the same, the strength of the compressibility 
shock can be varied (fig* 8). The locus of the thus obtainable 
states on the cone envelope, the connecting line of the end points 
of the individual meridian curves, is called "apple curve" because 
of its particular shape. Furthermore, the airspeed can be varied, 
so that to eaoh shock polar or each Mach nura.be.-. of flow there 
corresponds a different apple curve* 

The conical flold after the compressibility shock has a 
potential and is (u and o again denoting the velocity components 
in the axial section parallel and at right angles to the axis) 
characterized by the differential equation 


uu 


1 ♦ rf8 . ( 


U 


(7) 


The dashes indicate derivations with respect to uj or signifies 
the square of the local sonic velocity which for ideal gases is 
linked with the speed by the relation 

£ * (i£ * » 2 ) ' (8) 


The meridian curves are integral curves of equation (7) with the 
boundary condition (fig. 7) in the point ug, Pg of the shock polars 

with the slope of the connecting line from (uj, 0) toward (ug, Ug) . 
This means that the potential flow must precede the state of gas 
produced by the compressibility shock. Moreover, since the normal 
to the meridian curve indicates the polar angles of the respective 
state, the condition for the initial slope is necessary for reasons 
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of continuity. The end of the meridian line Is given hy the require- 
ment that the normal pass through the zero point of the hodograph: 
u + to 1 B 0. Direction of velocity and of normal must he identical. 

For the practical determination of the meridian ourve, a 
graphical method based on the shape (identical with equation (7)) of 
its different equation 


R ■ - 


JL 


-£ 


(9) 


is appropriate. Here (fig. 9) R is the radius of curvature, H 
the normal section between meridian curve ami U-axis, and U the 

u + op 1 


length of the perpendicular U - 

. Vi 

the hodograuh to the normal. As a 
(formula (8)), the curvature radius 


i- u 
is 
H 


<2 


from the zero point of 


a function of the velocity 
in a point of th e meridi an 

U, N, and w = \Ju + vr . 
progressive construction it 


in general, obtained fairly 


curve can be computed by marking off 
In the practical performance of this 
la found that the meridian lines arej 

accurately- This was checked by integration of equation (7) by the 
Runge-Kutta approximation process for ordinary differential equations- 
A departure in graphical accuracy occurs only at the meridian curves 
arising near the ends of the shock polars, that is near the axis. 

Most easily uncertain in the graphical determination is the position 
of the end point u^, ^ on the meridian curve, and with it the 

cone angle- In order to obtain it also relatively exact, individual 
meridian lines were mathematically unccked at points where their 
variation had been satisfactorily represented by the construction. 

The applicability even of the Runge-Kutta method to the meridian 
curves near the axis is limited- But here equation (7) affords 
approximations where the last term at the right Is taken equal to 

v* 

the constant initial value j> in the neighborhood of (uj_, 

Uj 2 ^ / a*2 ' \ 

and zrv in the vicinity of ( u-, , 0) . Quantities a. 

5 - . \ U 1 / 

aid are the sonic velocities related to (uj., 0) and (TI^, 0) 


0) 
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The shape of the meridian curves is then: 

*r.-lein^“ U 





supersonic 


(10a) 


*1 “ u / 

tt ■ tv. cos — — v flL - rrr? subsonic 

* « \T V 


(10b) 


The marked approximation character of these formulas is 
readily apparent frcrn the fact that In equation (10a) for all 
values of the parameter- A the curves pass through the point 


u ■ u-jy o-o, have the Bame Initial shape 


- 1 


and the 


Initial curvature aero. In equation (10b) the curves pass through . 
the initial points (u^, v>g) but have the same initial zero slope. 


In spite of that, the approximations are sufficient to reproduce 
the behavior of the apple curves near the "blossom" and the "stem," 
where a logarithmic singularity occurs: 


Ug « u^ 


*3 m \ 




"1 



subsonic 


(11a) 


( 11 %) 


On comparing the graphical Integration for the different shook polars 
it is found to be best for the maximum velocities. For the higher 
Mach numbers the curvature radii are too great to ensure good 
construction. However, this Is usually overcome by substituting. ' 
straight lines for the aros-. But for Mach numbers near unity, on 
the other hand, the construction becomes very inaocurate because - - 
of the abnormal curvature variation along the meridian line. For 
the vicinity of the critical velocity there is an approximation 
by Busemann which, however, does not become accurate except at 
very small shock polars as will be discussed later. 
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The airspeeds employed far the graphical numerical Integration 
of equations (7) and (9) correspond to the pressure factors (999.9915), 
9>6, 992, 986 , . . . 964, 956, 945, 930, 870.7 toy Busemana's oount 
In the Handbook far Experimental Physios (reference 5) ■ The subsequent 

n, 

table contains the Mach numbers M* ■ *3 or M ■ r— and the 

or a^ 

Mach angle referred to the critical and local sonic velocity for these 

°p 

pressure factors. The ratio of the specific heats it • — Is 1.405* 

°u 

The Mach numbers are a little more accurate than In the handbook* 

The bracketed pressure figure 999*9915 corresponds to a shock polar 
near tho critical velocity at which equation (7) was Integrated by 
the Busemann approximation method* 


TABUS I 


p- fact or 

U 1 

M* - -4 
a* 

u 1 

M « -=■ 
*L 

Mach angle 

dr 

(deg) 

Figure 

no. 

(999-9915) 

1.04106 

1.050 

79° 

13' 

*3 

* 

*- 

996 

1.1714 

1.2102 

55° 

10 ' 

33 

4 

10 

992 

1.2763 

1.3663 

47° 

2 ' 

46" 

5 

11 

988 

1*3662 

1.5045 

41° 

39' 

26 " 

5 

0 J 2 

984 

1*4472 

1.6403 

37° 

33' 

42" 

6 

13 

9B0 

1*5221 

1-7775 

• 34° 

14 1 

9" 

6 

14 

976 

1*5922 

1-91795 

31° 

25' 

31 

5 

15 

972 

1*6592 

2.0636 

28 0 

59' 

7 

7 

16 

968 

1*7206 

2.2157 

26° 

49' 


6 

17 

964 

1-7798 

2.3760 

24° 

53' 

22 " 

7 

18 

• 956 

1.8893 

2.7278 

21 ° 

30' 

20 " 

8 

19 

945 

2.0223 

3-3051 

17° 

36* 

42" 

8 

20 

930 

2.1729 

4-377 6 

: 13° 

12 ' 

19" 

8 

21 

©70-7 

2.4369 

00 

0 ° 
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Figures 10 to 22 represent the shock polars, meridian lines end 
apple curves for the several Mach numbers. (The number of the 
figure Is given in the last column of table I* ) The starting 
points for the meridian curves on the shock polars were so chosen 
that the plotted chords of the shock polars have constant difference 
of angle, A7 (second to the last column, table l) which amounts to 
4° to 9° depending upon the Mach number. For the related cone flows 
It Implies that at' the curve Mach number the half Included angles 
of the chosen .compressibility cone form a decreasing arithmetical 
sequenoe from 90° on. The lower limit is formed by the Mach number 
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which, however, does not. In general, belong to the sequence ■ The 
end jpolnto of the meridian- lines _ fix the. included angle of the cone 
in the' stream. They naturally form no simple sequence any more. 

The five lines plotted between shock polars and apple curve link 
the points of the individual, meridian curves in which "the difference 
between slope of normals and direction of speed is constant. With 
the notation of figure 7 they are the curves <p' - 1 * Constant. 

The value of the constants is noted at the individual lines, the 
numerals Indicate degrees, without it having been written in the 
figure. The direction of the normals (q>* ) can be computed far ' 
the points marked by the thin lines on a meridian line f.om the 
fairly accurately readable speed direction (^'). The reading 
of <p' from the figure is. In general, not quite so accurate* 

So for a good interpolation it is advisable to plot the relation- 
ship between f 1 and <p' in a right angled By stem of axis from 
the marked points, which then also gives, aside from the end points 
of the meridian line where <p* * V " Pi their starting point on 
the shock polars, where the direction o* the normals is at rlgjht 
angles to the plotted chord. For the meridian curves arising near 
the axis in the subsonic range the first apple curves carry only 
the points, the thin lines are emitted- The e:id points of the 
meridian lines further cany the ratio of the tank pressure before 
and after the compressibility shock, which naturally is already 
established by the point ug, u g on .the shock, polar. Moreover, 

the circle of the critical velocity + "i? «* a* ie indicated 

by dotted lines in each figure. 

The family of apple curves 1 b represented in figure 23. The 
lower part gives the shock polar for each apple- The inneimost 
dotted apple is obtained by approximation near the critical 
velocity a n d corresponds to tne bracketed digits in table I. The- 
thin system of curves permits a rough orientation of the points 
of the shook polars and’ apple curves given by the meridian curves- 
The exact ' relation is seen from figure 24(a), where the cone angle 23 
is plotted against the includod angle 27' of the compressibility 
cone for the individual shock polars- Figure 24(b) shows as 
contrast the corresponding relation in the two-dimensional case 
between wedge angle 25 and the angle of the compressibility 
shod; with respect to the flow direction 27. The aforementioned 
By stem of curves in figure 23 contains as first set the curves of 
equal entropy, that is, of fixed tank pressure ratio at the 
compressibility shock p 0 '/p 0 (formula (6)), as exemplified by 

■ 0-99, 0-9, and 0-7* The piece of the u-axie 1 « ^ * y * ■ j ^ 
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corresponde to the curve ■ 1... ,At. the border of the diagram, the 

p ° u i nr+i I> ' 

shock polar and apple curve for maximum velocity ~ * \ I — — ~i ■ 

a * V a “ 1 Po 

Ae second set the straight lines ug ■ Constant were chosen* By 
equation (5) the pressure riBe pg - p-]_ at shock referred to dynamic 
pressure jp^u^ 2 Is constant on the lines ug = Constant- The 

subsequent adiabatic pressure rise In the conical field after the 
shock up to pressure at the cone envelope is obtainable from 
the- gas states ug, vg and uo, o? by adiabatic compression of pg 
by the formula 


- Pi = P2 ~ Pi 

gPlUi 2 ^l u l 2 


It + 1 _ K - 1 + 

2 " 2 

R + 1 K - 1 U^ 2 
2 2 




P,. 


- 1 


Sk + 1 


1 ug^ + ^2 




a® 


"'I V 


* + 1 - !Ls - ,1 i. 

2 2 a*2 


Figure 25(a) shows the pressure ratio 


V3 

^iV 


plotted against cone 


angle for the several apple curves, that is, airspeeds* Figure 25(b) 


Po - 


represents the two-dimensional case of 


§ p i u i 


P 1 

Tf 1 against wedge angle* 


As in the 1 plane, there are, in general, two flows for every 
possible cone angle at fixed Mach number* On sufficiently small 
cone angles one solution Is subsonic, the other, supersonic, the 
latter occurring In practice* The range in which both are subsonic 
- ' is greater In this Instance than In the two-dimensional case* Here 
all points of the shock polar could be realized by considering a 
shock wave separated from the body* This is inapplicable here 
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'because there is so oonlcal field in the vioinity of an arbitrary 
-point of the shook wave* . Hence,, It remains- eut bpen question vhich 
range of the apple curveB and pressure . curves is of practical 
importance. The supersonic range should he assured. A study 
on slightly ourved oonlcal surfaces analogous to Crocco ’b investi- 
gation for the wedge is not yet available* 

The limit of the cone angle following from the apple curves 
for a specific Mach number is shown in figure 26, along with the 
maximum wedge angle* The angle sane 1 is characterized by the 
fact that for these angles wedge, as well as cone, flows exist, 
where the shock surface adheres to the edge ax tip. The second 
range includes the angles where the knife edges are worse, because 
the shook surface has become detached, but the tips are still useful* 
In the third range the shock surfaces for both the wedge and the 
oone lie upstream from the body. Figures 27(a) and 28(a) are 
modifications of figures 24(a) and 25(a)} the Mach number is 
consistently marked off on the axis of abscissa. The axis of the 
ordinate contains the size of the compressibility shock and pressure 
ratio. The individual ourvos correspond to cones of fixed included 
angles. The limiting curve 0 •* 0 shown as dash and dot line (-.-.-) 
in figure 26(a) indicates the pressure Jump before a blunt obstacle. 

If the adiabatic compression up to the stagnation at the body is 
also taken into account, it affords the plotted pitot pressure 
(Prandtl -Rayleigh) ■ 

Figures 27(b) and 28(b) represent the two-dimensional case. 

For several oone angles (supersonic range after the shock) the 
pressure calculation by Tlla) gives the simple approximation 
formula: 


% " Pi 
a p i“a 2 



at which the logarithmic singularity and Mach relationship is correctly 
defined, but the numerical factor 2 behind the logarithm is likely 
to be uncertain in certain circumstances* 


In conclusion, the integration of equation (7) for the apple 
curve in vicinity of the initial velocity is stressed. As demon- 
strated by Busemaxm, a suitable substitution makes it possible to 
range equation (7) into a differential equation Independent of the 
Mach number M and the specific ratio k. The seme substitution 
freoG the shock polar and hence the initial conditions of the 
meridian lines frcm. M and k. They are employed only 
for the determination of the end point of the meridian lines. 
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Figure 29 indicates the results of the integration in the coordinates 
x and 7 for the meridian lines Independent of M and it , which 
are related with u and o through ; - 

• ' Si 

y , S. +.A MLT. ft*. 

2 (M - l ) 3 / 2 a* 2 a* M *■ 1 

The initial abscissa on the shock polars Is written at the meridian 
curves’. The subsequent meridian lines (not shown) can be replaced" 

by curves 7 * ae tac . The related constants a. and b can be 
taken fr.am table II- 


TABLE II 


*2 

a. 

b 

-0.95 

0.00755* 

-3-59 

-•9 

.0298 

- 2.66 

-•5 

• 324- 

-1.48 

0 

•763 . 

-1-37 

• 25 

1.02 

-1.46 

■5 

1-39 

- 1-68 

•75 

2.28 

- 2.30 

• 9 

6.03 

-3-59 

• 97 

75-4 

- 6 . 6 5 


The end point of the meridian lines for a given M and k 
follows from the equation 

.21 x 

3 3 k + 1 (M - ir M - 1 


The dash denotes the derivation with respect to x. 

Translated b 7 J. Vanier 
National Advisory Committee 
for Aeronautics 
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Fig. 4a, b 



Figure 4a.- Schlieren photograph of wedge for adhering 
compressibility shock. 



Figure 4b.- Schlieren photograph of wedge for separated 
compressibility shock. 


NACA TM No. 1157 


Fig. 5a, b 







Figure 12. F-factor 988, 


Figure 13. P-factor 98/i 


U, u 

= 1,3662, — - 1,5015. 

a* a, 



Figure 14. 


P- factor 980, 


Hi 

a* 


1,5221, 


Figure 16 






si* 


Figur 


v 



Figur 
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Figs. 21,22 







NACA TM No. 1157 


Figs. 24b, 25a 



Figure 24b. Relationship between half wedge angle 8 and angle of slope between 
compressibility shock and flow velocity at different air speeds. 



Figure 25a.- Pressure on cone surface plotted against cone angle (3 

for different air speeds. (p« pressure at cone envelope, p 1 
pressure of flow, ^-P^u^ dynamic pressure.) 






Figs. 25b, 26 
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Figure 25b. Pressure jump pg - pj at compressibility shock referred to dynamic 

pressure in relation to wedge angle 6 for different air speeds. 



Figure 26. Limiting value of the half cone angle and wedge angle for adhering 
compressibility shock plotted against Mach number. 
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Figs. 27a-28b 



Figure 27a. Half included angle of Figure 27b. Slope TT of compressibility 

compressibility shock against the shock against flow for different wedge 

Mach number for different cone angles 2 S as function of the Mach - 

angles 2 ft . number. 



Figure 28a.- Pressure on cone 

surface in relation to Mach number 
for different included angles 2 p 


Figure 28b. Pressure jump at compression 
shock referred to dynamic pressure for 
several wedge angles 2 5 against the 
Mach number. 


(p Q pressure at cone envelope, p 1 

a 1 2 

pressure of flow , — p^u. dynamic 

pressure). 2 








